Dynamical control of two-level system's decay and long time freezing 
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We investigate with exact numerical calculation coherent control of a two-level quantum system's 
decay by subjecting the two-level system to many periodic ideal 2n phase modulation pulses. For 
three spectrum intensities (Gaussian, Lorentzian, and exponential), we find both suppression and 
acceleration of the decay of the two-level system, depending on difference between the spectrum 
peak position and the eigen frequency of the two-level system. Most interestingly, the decay of the 
two- level system freezes after many control pulses if the pulse delay is short. The decay freezing 
value is half of the decay in the first pulse delay. 
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I. INTRODUCTION 

It is of fundamental interests to study the quantum 
control of spontaneous emission which inspired the devel- 
opment of the quantum electrodynamics during the last 
century [H-Q ■ The unvifanted spontaneous emission often 
sets the ultimate limit of precise quantum measurement 
and many proposals have been made to suppress it 0-01 ■ 
One way widely applied to control the spontaneous emis- 
sion of a small quantum system, such as a neutral atom, 
is to place it either into or nearby a microcavity such 
that only a single mode or several modes of the cavity 
are resonant to the eigenfrequency of the quantum sys- 
tem [^] . In this way the structure of the vacuum is modi- 
fied and the spontaneous emission could be controlled by 
the properties of the cavity. Another way to control the 
spontaneous emission is dynamical control of the quan- 
tum system such that the coupling between the quan- 
tum system and the vacuum is effectively modified. For 
example, dynamical decoupling of the quantum system 
from the vacuum via either phase modulation or ampli- 
tude modulation could in principle extend the coherent 
time of the system [l^] ■ Analogous to the quantum 
Zeno effect (ZE) which says frequent measurements of a 
quantum system would prevent the decay of an unsta- 
ble quantum system |9l-[ll| , the extension of the coherent 
time through coherent modulations of the quantum sys- 
tem arc often called ZE as well. Under some unfavorable 
conditions, it is also possible that frequent modulations 
lead to acceleration of the spontaneous emission, the so- 
called quantum anti-Zeno effect (AZE) [l^ - fTsj . 

In practice people employ both ways separately or 
combination of them to realize the optimal control of 
a quantum system. G. S. Agarwal et al.^s work 0, [l^ 
is among many of such works. In their paper, a two- 
level model system with both structured vacuum and free 
space vacuum are investigated. They found significant 
suppression of spontaneous emission rate for structured 
vacua but either suppression or acceleration may appear 
in a free space vacuum, depending on the frequency of the 
control pulses which is 27r phase modulation pulses. They 



demonstrated for a free space vacuum that ZE shows up 
for wqt = 1 while AZE appears for ujqt = where ujq 
is the eigenfrequency of the two-level system and r is 
the delay time between control pulses. At the end they 
argued that AZE is possible for r > w^^^ 0. 

A puzzle arises if one accepts the AZE condition be- 
cause neither ZE or AZE appear if one leaves the quan- 
tum system alone in which case obviously r ^ '^o'^- 
this paper we revisit this problem by adopting the ex- 
act solution of a two-level quantum system which does 
not require the weak coupling and short time approxi- 
mations compared to Ref. Q- By investigating several 
typical structured vacuum, we show the conditions for 
quantum ZE and AZE and the boundary between them. 

Another puzzle is the ZE effect of a large number of 
pulses. According to the ZE, at a fixed evolution time 
t, the survival probability of the initially excited state 
approaches 1 with infinite number of pulses N ^ oo (the 
pulse delay t = t/N approaches 0). It is unclear what the 
road map looks like as N increases. Two ways might be 
possible. One way is that the rate of decay depends solely 
on the pulse delay r (independent on N) and approaches 
zero if T ^ |16|. Another way is that the decay rate 
depends only on the number of pulses N (independent on 
r) and the decay rate becomes zero after initial several 
pulses. For the latter one, we expect to see decay freezing 
at large number of pulses. We will find which way is the 
correct one in this paper. 

The paper is organized as follows. Section |TT] reviews 
briefly Rabi oscillation of a two-level system under a per- 
turbation and Sec. IIIII develops an exact formula of Rabi 
oscillation under periodic pulses. In Sec. IIVI after estab- 
lishing a quasi-level picture of the pulsed two-level sys- 
tem with constant spectrum intensity, we investigate the 
quantum Zeno and anti-Zeno effect for three spectrum 
intensities, including the Gaussian, Lorentzian, and ex- 
ponential one. The boundary between quantum Zeno 
and anti-Zeno effect is given. We discuss the decay freez- 
ing as pulse delay r getting small in Sec. El Finally, 
conclusion and discussion are given in Sec. I VII 
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FIG. 1: Diagram of a coupled two-level system subjected to 
ideal 2n pulses. 



II. RABI OSCILLATIONS OF A TWO-LEVEL 
SYSTEM 

Let us consider a simple model of a two-level quan- 
tum system with constant coupling to its environment 
modes. Such a model serves as the base of further discus- 
sions for more complicated systems with mode-dependent 
coupling. The Hamiltonian of the system with detuning 
6 = uj — ujo (we set fi = 1 for convenience) and coupling 
strength i/ is described by 0, [13] 

H = ^(|e)(e|-|5>(.g|) + ^(|e)(.9| + |5)(e|), (1) 

where ujq — — Eg is the difference between eigenfre- 
quencies of the excited state |e) and the ground state \g), 
with respectively eigenenergy Ee and Eg, and a; is the 
frequency of the external mode coupled to the two-level 
system [18'|. Utilizing Pauli matrices, the Hamiltonian 
can be rewritten as 

H = 5'a, + va^ (2) 

witha^ = \e){e\-\g){gla^ = |e)(.g|-f |5)(e|, and ,5' = 5/2. 
The evolution operator of the coupled two-level system 

is 

U = d,\e){e\+dl\g){g\+d2{\e){g\ + \g){e\) (3) 

where di — cos fit — i(5' /^l)smQ,t and ^2 = 
—i{v/^) sinili with il^ = 5'"^ + v"^. The transition prob- 
ability at time t from |e) to \g) is 

Peg{t) = ^Sin2(m) (4) 

providing the initial state is V'(O) = |e). The system re- 
vives at times such that ^It = ki: with k being an integer. 

III. CONTROLLED RABI OSCILLATIONS OF A 
TWO-LEVEL SYSTEM 

By applying an ideal 2tt pulse (or parity kick, see 
Fig. [ij , which is very strong in amplitude and short in 
time but gives a tt phase shift solely to the excited state 
|e) utilizing an auxiliary state \a) the state of the 
system changes according to 

x\e)+y\g) 27r pulse -x\e) + y\g) (5) 



where |a;p-j-|?/p = l. We denote such a pulse as Z pulse 
hereafter 

Z{x\e)+y\g)) = -x\e)+y\g). (6) 

In fact, Z ^ -cr^ = \g){g\ - |e)(e|. 

The evolution operator for the 2tt pulse and the free 
evolution t is 

u<" -zu^(X-k)' 

For N such operations, the evolution operator at time 
t = Nt becomes 

Let di — dir + idii and d2 — id2i. After some straightfor- 
ward simplifications with the use of Pauli matrices, one 
easily obtains 

U(N) ^ i^lcosNX-iaaSmNX) (9) 

where sinA = ~ ^^'-^ '^'^ ~ (1/sin A)(o?i,a2 — 
d2i(Jy). Note that A depends on r instead oit = Nt. 

The transition probability from the initial state |e) to 
the ground state \g) is 

Pe,{t = Nr) = |(5|C/W|e)|2=p,,(r)^^^.(10) 

sm A 

The above result is exact for any coupling strength and 
pulse delay r. For weak coupling v <^ 5, Q, ~ 5' and 
A ~ 5't + 7r/2 + 2k'K with k an integer, then 

Pea ~ TTTtan^^'rsin^ A^(5't (11) 

for even N , which is exactly Eq. (9) in Agarwal et a/.'s 
paper IS]. 

IV. QUANTUM ZENO AND ANTI-ZENO 
EFFECTS 

The above results are applicable only to single mode 
bath/environment which couples to the central two- level 
system. In general the bath has multimodes and the cou- 
pling may depend on the mode, e.g., the dipolar coupling 
between an atom and an electric and magnetic field. For 
a many-mode bath, the transition probability of the two- 
level system subjected to control pulses is in general given 

by 

with n the bath mode index [l^. Note that J7„ and A„ 
become mode dependent in the many-mode bath case 
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FIG. 2: Control pulse effect with fixed pulse delays, 
parameters are uo = 1, fo = 0.001, cJc = 100, and r = 
10, 4, 3, 2, 1 from top to bottom for solid lines. The dashed 
line denotes the free decay. Time is in units of l/uJo here- 
after. 



and V is constant here but will be considered as mode 
dependent in Sec. IIVDI Assuming the bath spectrum is 
dense, we turn the summation over mode index n into 
the integration over mode frequency cj, 



Peg = / dujp{uj)—siv?nT 



sin^ A^A 
sin^ A 



(13) 



where p(a;) is the density of states of the bath. For N = 1 
we get the free decay results (no control pulse). Before 
considering the real physical systems which usually have 
mode dependent coupling, we study several toy models 
with constant coupling strength v = v^y for all bath modes 
to gain some ideas about the effect of control pulses. 



A. Uniform spectrum intensity 

Taking p{uj) = po = I/lOc if o; e [0,a;c] and p{uj) — 
otherwise with lOc denoting the cutoff frequency of the 
bath, Fig.[2]shows typical free decay and controlled decay 
of the excited state. Except at very short times. Fig. [2] 
shows that the transition probability linearly depends on 
the total evolution time. Utilizing the linear dependence, 
one defines decay rate (Einstein constant) at long time t 



dt 



271-/90 t'o 



for free decay [?! and 

1 dp'^g _ 2itpqvI 



A' 
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i/2 + (Wfc - wo)V4 



(14) 



(15) 



for controlled decay, where ujk is the fcth resonant mode 
frequency [l^. For weak coupling vq ^ tt/t, the fcth 
resonant mode lies approximately at uJk ~ ujq + {2k + 
l)7r/r where sin A <C 1 and sin2AfA/sinA « TrS{X) w 
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FIG. 3: (Color online) Control pulse effect vs number of pulses 
(main panel) and relative positions of resonant modes to 
ojo and Uc (upper three panels A, B, and C corresponds to 
three regions in the main panel, respectively). The parame- 
ters are Nt = 100, = 1,vo = 0.001, and uJc = 100. Two 
big steps shows two major peaks out of the lower and upper 
cutoff frequency, respectively. These small steps shows minor 
peaks going out of the cutoff frequencies. Circles are for the 
frequency-dependent coupling case. 



{2h/t)5{uj — ojk) for large A^. Note that the two nearest 
neighbor peaks around wq with k = and k = —1 con- 
tribute equally about 40% among all the peaks 0, UBl ■ 

If we assume the bath has only positive frequency and 
concentrate on the two dominant peaks, we find 



~A 



0, UJcT < TT 
1/2, UqT < TT < 

1 , otherwise 



(16) 



and Peg/ Peg also shows step-like behavior as depicted 
in Fig. |3l Moreover, Fig. |3] exhibits only quantum 
Zeno effect, i.e., suppression of the decay by the control 
pulses. The decay are completely suppressed. A' ^ 0, 
once the control pulse frequency is larger than the cutoff 
frequency. 



B. Gaussian spectrum intensity 

A Gaussian spectrum intensity has a form as 
1 



p{u) 



12-kV 



/^r^, w>0 (17) 



where F and cj„i denotes the width and the position of 
the maximal intensity, respectively. 



1. OJo — 

In this case we expect only quantum Zeno effect to 
appear, A' < A, because it is easy to check that 

A 



7r2 ^ {2k -t- 1)2 p(c^o) 



< 1 



(18) 
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FIG. 4: (Color online) Reduced transition probability at fixed 
time (t — Nt = 100) vs number of pulses for Gaus- 
sian spectrum intensity for oJo = uJm (red dashed line) and 
ljJq = uJm — 2 (blue solid line). Circles and crosses are for the 
corresponding frequency-dependent coupling case. The ver- 
tical dot-dashed line shows the boundary between quantum 
Zeno and anti-Zeno effects for ljq — — 2. Other parameters 
are uo = 0.001, lJc = 100, cj™ = ajc/2, T = 1. 



where we have used the large TV and weak coupling as- 
sumptions. By considering the two dominant peaks, we 
further obtain that 



~A 



(19) 



for small r such that Fr < tt. Equation ((T^ shows that 
the decay rate decreases rapidly with the control pulse 
frequency in a Gaussian form. For exceedingly small 
T which satisfies Fr <C tt, A' is essentially zero, which 
means the transition is inhibited and the survival prob- 
ability of the initial state saturates. The red dashed line 
in Fig. m demonstrates the quantum Zeno effect and the 
prohibition of the transition. 



Strongest anti-Zeno effect is obtained at jwo — i^m\T = 
where one of the two main resonant modes lies near . 



C. Lorentzian and exponential spectrum intensity 

The spectrum intensity for Lorentzian and exponential 
are, respectively. 



p{uj) = 



F/tt 



(C^-CJ,„)2+F2 



w > 



(21) 
(22) 



Similar to Gaussian spectral intensity, only the quantum 
Zeno effect is observed if wo = w„i for the Lorentzian and 
exponential shape (Fig. [5]) . The transition probability is 
essentially inhibited once t <^ Tc where Tc ^ Once 
|wo — Wm| > F, we find both quantum Zeno and anti-Zeno 
effects and the boundary between Zeno and anti-Zeno 
regime is determined approximately by 2\ujm — ujo\t 
TT. The peak position of the quantum anti-Zeno effect 
lies about at l^o — ^^m\T ~ tt. As shown in Fig. [Sj the 
reduced transition probability in the exponential case has 
a narrower peak width than that in the Lorentzian case. 
Two small secondary peaks are noticeable in the anti- 
Zeno region of the exponential spectrum case. These 
peaks are due to the second and third quasi-level resonant 
to UJm- They satisfy the condition, respectively, \ujo — 
Wm|T « Stt, 57r. 



D. Frequency dependent coupling 



2. 



UJo 



> r 



For jwo — ujm\ < r, the results are similar to ujq = lo„^ 
and exhibit only quantum Zeno effect. To observe sub- 
stantially enhancement of the transition, i.e., the quan- 
tum anti-Zeno effect [l^, jwo — Wm| ^ T is required. More 
specifically, by taking the biggest peak around wqj one 
obtains the necessary condition for quantum anti-Zeno 
effect as 



|wo 



> 



41n - F w 1.3F. 



(20) 



The blue solid line in Fig. |4] shows the quantum Zeno 
and anti-Zeno effects for \ujq — Wml = 2. Both quantum 
Zeno (large A'' region) and anti-Zeno (small N region) 
effects are observed. At fixed time t = Nt, large r (small 
N) gives anti-Zeno effect while small r (large N) gives 
Zeno effect. The boundary between quantum Zeno and 
anti-Zeno effect is about 2\LUm—uJo\T ss tt if jwm— wq] ^ F. 



The widely adopted dipolar coupling in spontaneous 
emission of a two level atom or molecule has a frequency 
dependence as i^(w) = vq-Juj where vq is taken as a con- 
stant, which depends on the dipole matrix elements 
We will consider the same coupling spectrum intensity as 
before, i.e.. 



1 1 



w V27r F 
1 F/tt 



CJ (w - W„i)2 + F2 

1 1 



a; > (23) 
w > (24) 



a;2F 



(25) 



for Gaussian, Lorentzian, and exponential density of 
state, respectively. As shown in Figs. S] and [51 the fre- 
quency dependent coupling has little effect on the per- 
formance of the control pulses. 
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FIG. 5; (Color online) Same as Fig. 3] except that the spec- 
trum intensity is Lorentzian (top) and exponential (bottom). 

V. DECAY FREEZING 

By inspecting Eq. (1131) , the controlled transition prob- 
ability freezes if NX ^ 1 and the freezing value is 

1 f 1 

lim Peg ^ t; dojp{oj)—;Sm^nT—^, (26) 

where we have replaced the rapidly oscillating integrand 
sin^(A^A) with its average 1/2. The left column of Fig. [B] 
from numerical calculation indeed shows that the tran- 
sitions freeze at long times for three different cases at 
small T. The smaller the r is, the smaller the freezing 
value is. Freezing of the transition probability is equiva- 
lent to freezing of the survival probability of the excited 
state, f = l-p'^g. 

Moreover, for small enough pulse delay r, we have 
sin^ Awl thus 

1 f 1 

The above relation indicates that the freezing value of the 
transition probability at long times is one half of the free 
transition value at i = r which is the pulse delay. Clearly, 
the road map to the quantum ZE is that the decay rate 
becomes zero after initial pulses (decay freezing) and the 
freezing value of the survival probability / approaches 1 
as T decreases. 

Define differential freezing value sir) = "p'^giNT — 
1000) — \peg{T), which describes the difference of the 
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FIG. 6: Decay freezing (left column) and dependence of the 
differential freezing value on pulse delay (right column) for the 
Gaussian (top row), Lorentzian (middle row), and exponential 
(bottom row) spectrum. Only results of odd number pulses 
are plotted for a better view. The dashed lines show decay 
without control pulses. Crosses of dashed lines and solid lines 
in the left column panels denote the specific pulse delay r. 

transition probability after many pulses and one half of 
the decay in the first pulse delay. The right column 
of Fig. ini shows that e decreases with increasing of 1 /r 
(decreasing of t), confirming the analytical results of 
Eq. (EZD. 

In fact, decay freezing exists not only in the model 
system we consider in the paper, but also exists in many 
other pulse-controlled systems, such as gated semicon- 
ductor quantum dotpO^lll , spin-boson model [1^ [1^ , 
and nuclear spins [23|. The basic idea behind the de- 
cay freezing is that the control pulses creates an effective 
preferred direction along which the decay is frozen. In 
terms of Pauli matrix (c.f. Eq.[2]), the preferred direction 
created by control pulses in the model we consider is z. 

VI. CONCLUSION AND DISCUSSION 

By investigating three models of coupling spectrum 
intensity (Gaussian, Lorentzian and exponential), we 
demonstrate that a two-level system subjected to many 
ideal 27r pulses exhibits both quantum Zeno and anti- 
Zeno effect, depending on the relative position of loq to 
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the peak position Wm of the spectrum and the pulse delay 
r. Instead of decreasing the decay rate, the pulsed two- 
level system shows decay freezing after many pulses at 
small r and the freezing value of the survival probability 
of the initial excited state approaches 1 (no decay) with 
decreasing r. 

In this paper, all the spectrums have single peak and 
we observe only single quantum Zeno and/or anti-Zeno 
region. Under some special circumstances where a mul- 
tiple peaks spectrum exists, one would expect multiple 
quantum Zeno and anti-Zeno regions. We have also as- 
sumed that the spectrum of the structured vacuum is 
time independent where the back action exerted on the 
vacuum by the two-level system has been neglected. A 
full quantum version of the coupling between the two- 
level system and the vacuum could possibly change the 
picture of the controlled decay at long times but the 
short time behavior would be intact in the weak cou- 
pling regime, because the back action is weak and needs 
a long time to manifest its effect on the two-level system 
dynamics. 

We consider only periodic pulse sequence (r is fixed) 



in this paper. In principle, other pulse sequences with 
varying r may also give similar quantum Zeno and anti- 
Zeno effect. They may even have additional advan- 
tages [2^ [2^. In addition, the strength and the dura- 
tion of the 27r pulses are finite in practice. One could 
minimize the finite pulse effect by employing the phase 
alternation techniques [2^ or the Eulerian protocols [13 ■ 
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